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As a prototype model of antiferromagnetism, we propose a repulsive Hubbard Hamil- 
tonian denned on a graph A = AVJB with AnB = and bonds connecting any element 
of A with all the elements of B. Since all the hopping matrix elements associated with 
each bond are equal, the model is invariant under an arbitrary permutation of the 
,4-sites and/or of the ,6-sites. This is the Hubbard model defined on the so called 
(Na, A^b)- complete- bipartite graph, Na (Nb) being the number of elements in A (B). 
In this paper we analytically find the exact ground state for Na = Nb = N at half 
filling for any JV; the repulsion has a maximum at a critical iV-dependent value of the 
on-site Hubbard U. The wave function and the energy of the unique, singlet ground 
state assume a particularly elegant form for N — > oo. We also calculate the spin- 
spin correlation function and show that the ground state exhibits an antiferromagnetic 
order for any non-zero U even in the thermodynamic limit. We are aware of no previ- 
ous explicit analytic example of an antiferromagnetic ground state in a Hubbard-like 
model of itinerant electrons. The kinetic term induces non-trivial correlations among 
the particles and an antiparallel spin configuration in the two sublattices comes to be 
energetically favoured at zero Temperature. On the other hand, if the thermodynamic 
limit is taken and then zero Temperature is approached, a paramagnetic behavior re- 
sults. The thermodynamic limit does not commute with the zero- Temperature limit, 
and this fact can be made explicit by the analytic solutions. 



I. INTRODUCTION 



The Hubbard Hamiltonian is one of the most popular models to describe strongly correlated electron 
systems. In spite of its simple definition, it can be exactly solved only in few cases. An example is the 
milestone solution of the Hubbard ring [1] by means of the Bethe-Ansatz [2] extended to fermions [3] [4]. 
Other exact solutions were worked out for the one-dimensional Hubbard ring in the presence of an external 
magnetic field [5], for the one-dimensional Hubbard chain with open boundary conditions [6] and thereafter 
for the SU(N) one-dimensional Hubbard ring [7]. However, when the space dimensionality becomes bigger 
than 1, few exact results are available and usually they concern the ground state properties. Among them, 
we mention the Lieb theorem on the ground state spin degeneracy [8] that will be explicitly used in this 
work. Exact ground state wave functions are even more infrequent. To the best of our knowledge the only 
non-trivial results are the ferromagnetic ground-state solutions devised by Nagaoka [9] (infinite repulsion and 
one hole over the half filled system), by Mielke and Tasaki [10] [11] (for which the kinetic energy spectrum is 
macroscopically degenerate or at least very flat) and by Wang [12] (infinite-range hopping). 

In the Hubbard Hamiltonian with infinite-range hopping a particle in any site can hop to any other site of 
the system; the associated graph is said to be complete. This model was numerically studied by Patterson on 
small clusters [13] in 1972 and solved in the thermodynamic limit by van Dongen and Vollhardt [14] only at the 
end of the eighties. Much more effort was needed to find the exact ground state(s) in the finite-size system. 
Verges et al. managed to accomplish this task for arbitrary numbers of particles and sites in the limit of 
infinite on-site repulsion U [15] by exploiting a scheme proposed by Brandt and Giesekus [16]. Two years later 
Wang constructed explicitly the ground states of the system for arbitrary U and number of particles above 
half filling [12]. In the case of one particle added over the half-filled system the ferromagnetic ground-state 
solution follows as a special (and the easiest) case of general results by Mielke and Tasaki [17]. 

In this paper we find the exact ground state wave function of the half-filled Hubbard model on the Complete- 
Bipartite Graph (CBG) for arbitrary but finite U. The CBG A = A U B has bonds connecting any element of 
A with all the elements of B and can be considered as the natural further step (with respect to the complete 
graph previously described) towards the standard Hubbard Hamiltonian (defined on the hypercubic lattice and 
hopping between nearest-neighbours sites). Even if the CBG is still somewhat unrealistic, we feel that exact 
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solutions should always be welcome, expecially because they lend themselves to be generalized. Furthermore, 
our solution is an example of antiferromagnetic ground state in a model of itinerant electrons, in contrast with 
the ferromagnetic solutions mentioned above; it may provide useful hints about antiferromagnetism outside 
of the strong coupling regime (where the Hubbard model can be mapped onto the Heisenberg model). 

The paper has been organized as follows. In Section II we define the model and discuss the physics of 
the non-interacting (U = 0) Hamiltonian together with few relevant examples of finite-size realizations. In 
Section III we study the thermodynamic limit of a class of Hubbard-like models, which includes the Hubbard 
Hamiltonian on the complete graph and on the complete bipartite graph, having a non-extensive number of 
isolated one-particle energy levels plus a single level whose degeneracy is proportional to the size of the system. 
Following the reasoning of van Dongen and Vollhardt [14], we show that the kinetic term is totally decoupled 
from the interaction term and that the system behaves as a paramagnet: the spin-spin correlation length is 
zero. These results are a consequence of the trivial behaviour of such models any time the thermodynamic limit 
is taken first. On the contrary, much more difficult is to find exact properties in finite-size systems and different 
answers may be obtained if the thermodynamic limit is taken only at the end of the calculations (as we shall 
show in Section VI). In Section IV we introduce the essential tools to face the problem. Let N = Na = Nb 
be the number of sites of each sublattice A and B. First, we find a (2JV — 2)-body determinantal eigenstate 
|$af) of the Hamiltonian with vanishing double occupation; then, we demonstrate that it is a key tool to 
build the ground state. We show that mapping the .4-sites onto the S-sites and viceversa, \$af) retains its 
form except for a spin-flip; we shall call this property the antiferromagnetic property. Several analogies with 
the properties of the half-filled standard N x N Hubbard model are pointed out at this stage. Here we also 
deal with the spin projection of onto the singlet and the triplet subspace and useful identities between 

the two spin-projected states are obtained. Then, in Section V we propose an Ansatz for the ground state 
wave function at half filling containing the singlet and triplet projections of \<&af)- We set up the Schrodinger 
equation and by exploiting the antiferromagnetic property we close the equations and get 3 exact eigenstates. 
The ground-state uniqueness proved by Lieb [8] is used to show that the lowest-energy state of our Ansatz 
actually corresponds to the ground state of the system. Remarkably, the ground state energy is negative 
for any value of the repulsion U ; qualitatively, we may say that the particles manage to avoid the double 
occupation very effectively. In Section VI we study the ground state energy Eo as a function of U and of the 
volume of the system and we discuss the implications of exchanging the thermodynamic limit with the limit 
of zero temperture. We find that Eo is a monotonically increasing function of U and N due to the existence 
of non-trivial correlations even for large N. The nature of these correlations is investigated by computing 
the expectation value of the repulsion. We show that for any finite iV there is a critical value of U yielding 
maximum repulsion. From the exact spin-spin correlation function we find that the ground state average of 
the staggered magnetization operator squared is exstensive. Hence, the ground state is antiferromagnetically 
ordered; we underline that this holds not only at strong coupling, but for any value of U. Finally, a summary 
of the main results and conclusions are drawn in Section VII. 



II. DEFINITION OF THE MODEL 



Let A = A U B with A H B = be a collection of sites and 

\A\=N A , \B\=N B ^\A\=N A + N B . (1) 

Here and in the following we shall denote by |«S| the number of elements in the set S. We consider the Hubbard 
Hamiltonian 

-^Hubbard = H + W, H = ^ y]tx, V ct a Cy<r, W = lj'^h x} h xl (2) 

where c xa (c|, CT ) is the annihilation (creation) operator of a particle at site x with spin a =f , { and n xa — c xa c x(y 
is the corresponding particle number operator. The hopping matrix T with elements T XtV = t XtV = t y , x is a 
real-symmetric matrix while U is a positive constant. If 

t for x G A (x G B) and y G B (y G A) , , 

otherwise 

the graph A is said to be complete bipartite and we call the Hamiltonian in Eq.(2) the CBG-Hubbard Hamil- 
tonian. The model is invariant under an arbitrary permutation of the ,4-sites and/or of the £>-sites. Therefore, 
the symmetry group includes Sn a <8» Sn b , Sn being the set of permutations of N objects. As usual, the full 
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group is much bigger: the presence of spin and pseudospin symmetries [8] [18] leads to an SO(4) internal 
symmetry Group [19] [20] and in the case Na = Nb there is a Z2 symmetry because of the A <-> B exchange. 

In Fig.(l) we have drawn a few examples of finite-size systems. For N = 1 and N — 2 the model is equivalent 
to a one dimensional ring of length L — 2, 4 respectively. For N — 3 we have what can be understood as a 
prototype, (1,1) nanotube model, the one of smallest length L = 1, with periodic boundary conditions. For 
general N, one can conceive a gedanken device, like the one illustrated pictorially for N = 6 in Fig.(l.d). 
(The whole device should actually have the topology of a torus, and the two horizontal faces should coincide, 
but this is not shown in the Figure for the sake of clarity.) The N vertical lines represent a realization of 
the A sublattice while the B sublattice is mimicked by the central object. The radial tracks in the Figure 
represent conducting paths linking each A site to each B site according to the topology of our model. The A 
sites are represented by one-dimensional electron-in-a-box systems of length L. Each A site has the one-body 
energy spectrum Ef roc ~ n 2 /mL 2 where n is an integer and m is the electron mass. We assume that L is so 
small that the excited states are at high energy and can be disregarded in the low-energy sector; this requires 
ml 2 < l/U. Here, of course, U denotes the Coulomb self-energy of a box with two electrons. The B sites 
are hosted by N quantum dots, that can be represented by <5-function-like attractive potential wells of depth 
V; if \V\ is large and the radius of the dots is <C L, these are practically independent of each other. The 
one-body energy of each B site is Ebound ~ —mV 2 ; the unbound states can be neglected if mV 2 3> U. We 
are assuming for simplicity that the U of the B sites is the same as for the A ones. Turning on a constant 
potential Vo on the central object, one can arrange that the Hubbard Hamiltonian on the CBG -Hrrubbard gives 
a good description of the system, with filling < 2 per site. 



t N=6 




FIG. 1. Physical representation of the CBG for N = 1 (a) and N = 2 (b) which is equivalent to a one dimensional 
ring of length L = 2, 4 respectively. For N = 3 (c) we have the (1,1) nanotube of smallest length and periodic boundary 
conditions. For TV = 6 (d) we present the gedanken device described in the text. 



It is convenient to label each site with an integer in such a way that x = 1, . . . , Na corresponds to sites in 
A and x = Na + 1, • • • , Na + Nb corresponds to sites in B. In the special case Na = Nb = N the hopping 
matrix can be written as 



Z J 
J z 



(4) 



where Z is the N x N null matrix and J is the N x N matrix whose generic element J X:V = 1. The one-body 
spectrum has three different eigenvalues 

e g = ~Ni= -t 
so = 

Eg =Nt = t (5) 

with degeneracy d g = 1, do = 2N — 2 and dg = 1 respectively. We use the convention t > so that e g is the 
lowest level and we shall call Shf the set of zero-energy one-body eigenstates. 
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The orthogonal matrix O that diagonalizes T can be written in the form 



O = 



V 



i i 

Pin V2N 
R 

Zy 
1 1 

/2JV V2JV 



-1 -1 

/2N V2N 

R 

_L 1 

/2N V2N 



(6) 



/ 



where Z r is the (N —1) x N rectangular null matrix, while R is an (N — 1) x N rectangular matrix whose rows 
are orthonormal vectors which are orthogonal to the A— dimensional vector (1, 1, . . . , 1, 1). The zero-energy 
one-body orbitals can be visualized by a simple argument. Consider any pair x, y, with x ^ y, of sites 
belonging to the same sublattice, say A, and a wavefunction ip x , y taking the values 1 and -1 on the pair and 
elsewhere in A and in B. It is evident that tp x ,y belongs to Shf- Operating on ip x ,y by the permutations of 
Sn we can generate a (non-orthogonal) basis of N — 1 eigenfunctions [21] vanishing in B; further, by means of 
the Z2 symmetry, we obtain the remaining orbitals of Shf, which vanish on A. This exercise shows that the 
group considered above justifies the (2N — 2)-fold degeneracy of the one-body spectrum. 

Having obtained the one-body basis, one can form a suitable R matrix by orthogonalization. Let us define 
the transformed operators: 



9 = ^ Oi, x c x = ^_ ^(cz — c x+N ) 

x = 1 x — 1 

2JV JV 

ai = Oi+i, x Cx = ^ Ri, x c x , i = l,... ,N — 1 

a;— 1 x — 1 

2JV JV 

h = Oi+N,xC x — Ri,xC x +N, i = 1, ■ ■ ■ , N — 1 

a;— 1 x — 1 

2JV JV 

3 = J]] C^JV.xCa; = y^(Ca; + Ce + jv). 



(7) 



The inverse transformation reads 



Cx = -^7=(g + 9) + ^ Rj,x<n, xeA 

i 

c x = -y=(-g + g) + y]Ri,x-Nbj, xeB 



and the kinetic term Ho becomes 



Ho 



(8) 



(9) 



Hence, if we do not rescale the hopping constant, the average kinetic energy remains an extensive quantity 
proportional to Nt: the kinetic energy of the two particles in the lowest level coincides with the kinetic energy 
of the whole system. On the other hand, if t ~ 1/N the two energy gaps £0 — £ g and Eg — so remain constant 
as N increases. As we shall show in the next Section, the model can be exactly solved in the thermodynamic 
limit whatever is the dependence of e g and e s on |A|. 



III. THERMODYNAMIC LIMIT 



The Hubbard Hamiltonian on the CBG (as on the complete graph) belongs to a class of Hubbard-like 
models that can be exactly solved in the thermodynamic limit. Here we generalize the results by van Dongen 
and Vollhardt [14] to the case of a non-extensive number of non-vanishing one-body eigenenergies and in the 
presence of a local external magnetic field coupled to the spin of the particles. Even if the reasoning is similar 
to that of Ref. [14], a detailed presentation of the main results is needed to clarify more subtle points and 
make our work self-contained. Furthermore, some of the results are absent in the original paper and we believe 
that their derivation will facilitate the reader in the comprehension of what follows. 
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Let us consider the Hamiltonian 

H — //Hubbard — fJ- ^ ] ^ ' fixer — ^ ' /^(^rcT — n x \) = -Hliubbard + -f/p + Hh (10) 

where -^Hubbard = -f/o + W is defined in Eq.(2), \i is the chemical potential and H h represents the coupling 
with an external magnetic field. Let us assume that only a finite number f of eigenvalues of the hopping matrix 
T are non-vanishing, e» =fc for i = 1, . . . ,f, while the remaining A| — f arc identically zero. Denoting with <j>i 
and <j>\ the annihilation and creation operators that diagonalize Ho we have 

{ 

i — 1 <7 

and the CBG corresponds to the case f=2, see Eq.(9). The Gran-Partition Function Zp [T, U, {h x }] can always 
be written as 

Z fj [T,U,{h x }} = Zp[T,0,{h x = }}e w ^ T > u >^\ Wp[T, U, {h x }} = (e~ So *'M" T )g (i 2 ) 

where j3 is the inverse Temperature, Hi(r) = W(r) + Hh{t) is the interaction Hamiltonian in the imaginary- 
time Heisenberg representation and (...)§ is the thermal average in terms of Zp [T, 0,{h x = 0}] = Z° p [T] where, 
by the linked-cluster theorem, one has to retain only those contributions represented by connected diagrams. 
As far as f remains fixed, we can safely substitute W/3[T, U, {h x }] with W/3[T = 0, U, {h x }] in Eq.(12). Indeed, 

Wp[T,U,{h x }]=Wp[T = 0,U,{h x }] + O(i/\A\) (13) 

as long as the unperturbed one-body Green's functions G; ct (t) = (T0 iCT (r)0j CT (O))o, i = 1, ■ ■ ■ ,f, do not diverge 
in the thermodynamic limit; here T is the Wick time-ordering operator and (. . .)o is the thermal average with 
U — and {h x = 0}. From the explicit expression of G;<t(t), 



■< e «-"> T (l- fi) r>0 
-e-^-^fi t < 



with /j = [1 + e l3 ^ Si ~^]~ 1 the Fermi distribution function, it follows that Gi^ir) converges to a finite value 
whenever lim| A |^oo£i is we U defined. One can use the result in Eq.(13) to write the exact Gran-Partition 
Function in the thermodynamic limit 

Zp[T, U, {h x }] = z f°f} Q] MT = 0, U, {h x }] = 



=n 1+2ze n+^ e """ n[ i+2 — h ^+^ e ~H - ( i5 ) 

i=i L -I ieA 

where z = e 13 ^ is the fugacity, and the thermodynamic potential £lp[T, J7, {fea;}] = — lnZp[T, U, {h x }] 

Qp[T, U, {h x }] = -j- I 2f ln(l + z) - ln [! + 2^e _/3£l + zV 2/3ei ] - ^ ln t 1 + 2z cosh /3ft* + z V 317 ] J 



i6A 



(16) 



Eq.(16) reduces to the result of van Dongen and Vollhardt in the case f=l, h x = Vcc € A and ei oc — |A| < 0. 
All the thermodynamic quantities can be derived from Eq.(16) and we defer the interested reader to the 
original paper by van Dongen and Vollhardt for further details. Here we calculate the magnetization at site 
x, m x — — ^{dQ,p/dh x ) = \{n x <\ — h x {) = (S x ) and the connected spin-spin correlation function {S x Sy) c = 
(S x Sy) — (S x )(Sy) = —(l/4:[3)(d 2 Qf3/dh x dhy), where (. . .) means the average in the gran-canonical ensemble. 
Due to the exact decoupling of the kinetic energy, we expect from the begining a paramagnetic behaviour. 
We have 

1 dtlp _ zsinh j3h x 

mx ~~2 'dhZ ~ 1 + 2z cosh f3h x + z 2 e-l 3u [ ' 

where z can be expressed in terms of the number density by means of the relation 

1 dQp _ 2z ^ cosh f3h x + ze'^ 



|A| dn |A| ^ 1 + 2z cosh /3h x + z 2 e-P u 



(18) 
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Since z is finite for any non-zero Temperature, m x — > when h x — > 0, i.e. 



lim lim lim m x = . (19) 

/3^oo /i^^O |A|^oc 

On the other hand, m x / if we exchange the last two limits in Eq.(19). Indeed, taking h x = ft=const> 
for the sake of clarity, Eq.(18) can be exactly solved for z: 

e^ u + h )( n -l)/(2-n) n> 1 

e' 3C7 / 2 [n/(2-n)] 1 / 2 n = l . 

e~ 0h n/(l -n) n < 1 (20) 



_ yj(n - l) 2 cosh 2 /3ft + n(2 - n^-^ - (1 - n) cosh j3h 
Z ~ {2-n)e-' 3u 



Substitution of these asymptotic behaviours in Eq.(17) yields 

r (2-n)/2 n > 1 

lim lim lim m x — < 1/2 n = 1 . (21) 

0+fl-oo|AHoo y n/2 n<1 

The case ft < is similar and one can show that Eq.(21) remains true in the limit h x — > CP if m x — > —m x . 
Therefore, the magnetization does not depend on U and the zero- Temperature limit of the model is the same 
of the paramagnetic Hamiltonian i/ par a = + H h , where and H h are defined in Eq.(10). From Eq.(17) 
we also conclude that 

(sisir = (aSfc) = ^ {slsZy) = {s ' ){s5) = mxmy ' (22) 

a result which is independent of the local external field configuration: in the thermodinamic limit two localized 
spins do not interact and the spin-spin correlation length is zero. Denoting with G sp i n (a;, y; /3, ft, |A|) the thermal 
average of S x Sy at inverse Temperature /3, external field ft and size of the system |A|, from Eq.(19) we get 

lim lim lim (S x Sy) = lim lim lim G sp in{x, y; j3, ft, |A|) = , (23) 

/3^oo h^O |A|^oo (3— >oo h^O |A|— >oo 

that is, the spin-spin correlation function vanishes if we first take the limit ft — ► and then the limit of /3 — > oo. 
We emphasize that the above results are correct only if we first take the limit |A| — > oo. The thermodynamic 
limit makes the model trivial and different graph structures, like the complete graph and the complete bipartite 
graph, all yield the same thermodynamic behaviour. On the contrary, a much more difficult task consists in 
finding exact results in finite-size systems. To face such problems is not only a mathematical exercise. After 
that in Section V the exact ground state of the half-filled Hubbard model on the CBG is explicitly written 
down for arbitrary value of U and |A| and in the absence of an external field ft, we shall see in Section VI that 
the thermodynamic limit, |A| — » oo, and the limits of zero Temperature and external field, f3 — > oo, ft — > 0, 
do not commute. 



IV. THE W = STATES 



In this Section we introduce the essential tools to deal with the antiferromagnetic ground-state solution. Let 
us consider the one-body eigenstate a]\0) with vanishing amplitudes on the B sublattice, that is (0|ca;aj|0} = 
if x G B. Similarly, b\\Q) has vanishing amplitude on the A sublattice and therefore the (2N — 2)-body state 

|^>=aL...a t _ lCT 6^...6 t _ 1 ^|0>, a = -a (24) 

is an eigenstate of Ho and of W with vanishing eigenvalue. In the following we shall use the wording W = 
state to denote any eigenstate of //Hubbard in the kernel of W . It is worth to observe that by mapping the 
,4-sites onto the Z?-sites and viceversa, |$af) retains its form except for a spin- flip; we call this property the 
antiferromagnetic property for obvious reasons. 

In the non-interacting (U = 0) half filled system, the structure of the ground state is trivial: two particles 
of opposite spin sit in the lowest energy level g, while \Shf \ = (2N — 2) particles lie on the shell Shf of zero 

/ 2N — 2 \ 2 

energy. In the spin S z = subspace this ground state is I ^ J times degenerate. To first order in W, 

the degeneracy is only partially removed [22]. Indeed, if Ps,m s is the projection operator onto the subspace 
of total spin S with ^-component Ms, the structure of the determinantal state in Eq.(24) yields 

Ps,o|*&> = |*^>^0, VS = 0,...,JV-1; (25) 
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therefore, the states {gjgjl'l'^}, S = 0, . . . , N — 1} belong to the ground state multiplet in first or- 
der perturbation theory (being W = states). Since the Lieb's theorem [8] ensures that the interacting 
ground state \^q(N,U)) is a singlet, only can have a non- vanishing overlap with \^o(N,U)) in 

the limit U — > + . In Section V we shall prove that <?{g{| ( I ) ^ F ) is the ground state for U = + , that is 
\imu^o{^o(N, £/)|(/j<7j| ( I>^) = 1. In this model, the first-order solution has a peculiar significance. Usu- 
ally, the exact ground state has no overlap with the non-interacting one, in the thermodynamic limit: this 
is because the U = ground states define a proper subspace of the full Hilbert space. Below, we prove that 
{9q(N, U)\g^ g^l&Ap) 7^ ^U; that is, the lowest approximation keeps a finite weight. 

At this stage, we note the analogy of the above results with those relevant to the standard Hubbard model 
defined on a N x N square lattice with periodic boundary conditions (hopping matrix elements only between 
nearest neighbor sites). The determinantal state \<&^f) resembles the W = state that we obtained for even 
N in Ref. [23] [24] . The degeneracy of the zero-energy one-body eigenspace was 2JV — 2 and it was shown [25] 
that N — 1 zero-energy eigenfunctions have vanishing amplitudes on a sublattice while the remaining N — 1 
vanish on the other. 

As we shall see in the next Section, we need the explicit expression for the singlet and the triplet 

|'I > ^™), m = 0, ±1, to calculate the ground state of the CBG-Hubbard model. Therefore, we now briefly review 
how to get the projection of the determinantal state \&af) °nto the singlet and the triplet spin subspaces. 

To obtain the singlet |$^) one has to antisymmetrize each product ot^ftt^, getting a two-body spin singlet 
operator, and subsequently antisymmetrize with respect to the N — 1 indices of the 6^'s; one sees easily that 
this entails the antisymmetrization of the af's. Hence 

N-l 

i«=EH p n-U)i°> ^ 

p i=i 

where <j\ ■ creates a two-body singlet state: 

and P is a permutation of the indices 1, . . . , N — 1. Since W commutes with the total spin operators and the 
determinant \^af) is a W = state, = as already noted. Applying W on the state in Eq.(26), one 

can verify this property by direct inspection. 

In a similar way one obtains the triplet projection. We define y as the two-body triplet creation operator 
with vanishing ^-component m = in the orbitals ai and bj : 

= ^4Ai+4A)- ( 28 ) 

Then, one has 

N-l 

i^) = E(-) p E^4n-U)i°)- ( 29 ) 

P 3 = 1 ift 

The m = ±1 components of the triplet in Eq.(29) can be otained by means of the total raising and lowering 
spin operators. Let us introduce the following notations 

JV-l 

S+ a = al r a iU Sp = [S+ a ]\ S +a = J2 S ? a ' S~ a = [S +a ]* (30) 

N-l 

St b = b^b iU Si b = [S+ b ]\ S-+ 6 = ]T ^ S- b = [S +b ? (31) 

i=l 

N-l 

t r< za ^ / ^ a ^ a \ £. a \ y ^ a ^x a £. a i £. a o za ^ I £. a £. a \ /on\ 

n i* = aLa irT , S t =2( n iT _n ii)> n ° = n ^-> n = n T+ n i> S =^ T ^ ^ > 

i=i 

N-l 

nl=blb ia , Sf = \{n b ] -h b l ), n^ = E«t, n b = h^+h b , S* b = \{n\ - n\) (33) 

(34) 
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The states l*^ 1 ) can be written as 



P 3 = 1 i& 

l^- 1 ) = -^(sr- + fi r»)|*£> = E(-) P E T ^w t )II ff Uol > ( 36 ) 



3 = 1 i& 



with 



Equivalently, the triplet state |$^™ ), m = 0, ±1, can also be expressed in terms of the singlet |$^}. It is a 
simple exercise to prove the following identities 

N-l 

- S+ 6 )|<0 = (S +a - S +b )\<f°/ F ) = -V2\^/ F ) (38) 

JV-l N-l 

E(™a ~4)lO = E(-"U4)lO = {-S' a + S' b )\*°£) = -|<^°> (39) 

i=l i=l 
JV-l 

E(-*r + sr b )\*°j° F ) = (-s- a + s- ft )i<o = -^i^ 1 ) . (40) 

i=l 

As for the singlet, one can verify that W\$>\™) = 0, m = 0, ±1, using the definitions in Eqs. (29-35-36). 



V. THE GROUND STATE AT HALF FILLING 

We are now ready to calculate the ground state |^'o(^, U)) of the half filled CBG-Hubbard model described 
in Section II. As already observed, <?{<?{ | j?) is a good candidate for the non-interacting ground state; its 
quantum numbers are the same as those of |*o(7V, U)) and moreover it belongs to the first-order ground state 
multiplet. Using the short-hand notations 

|«; >^{«;t|0>, tf)=g\g\\0), \[gg}°> ) = ±={g\g\ - g\g\)\Q) (41) 

for the three different two-body singlets that one gets from the lowest and the highest energy orbitals g and 
g and 

Hgg] 1 ' 1 )^^), Hgg] 1 ' )^ -±=(g\gl+glg\)\0), fog] 1 -- 1 ) = g\g\\0) (42) 

for the triplet, we propose the following Ansatz for the interacting ground state |*o(A r , U)): 

1 

\MN,U)) = [j g \g )+j s \g")] ®|<0+7o E (-) m |[^] 1,m >®l*AF m > , (43) 

m= — l 

where the 7's are c-numbers. It is worth to note that in |\E , o(A r , U)) the number of particles in the shell Shf is 

a constant given by 27V — 2. A priori, there is no reason for this choice since the total number operator n a + n b 

of particles in the shell Shf does not commute with the Hamiltonian. Nevertheless, we shall see that the 

scatterings which do not preserve this number cancel out provided in S h / the (27V — 2)-body state is a W = 

state. We shall refer to this very remarkable property as to the Shell Population Rigidity. We emphasize that 

even constraining the particle- number in Shf to be 27V — 2 with vanishing total spin ^-component, there are 
\ 2 
27V - 2 1 

configurations that can contribute to the ground state expansion in the interacting case, while 



N-l / 

the Ansatz (43) contains only |$^) and |$^™), m = 0, ±1. The reason why the W = states \^^ F S ) with 
S > 1 do not enter into the ground state expansion (43) comes from the Lieb's theorem [8]: the ground state 
must be a singlet and with only two particles outside Shf (in the g and/or g states) the angular momentum 
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composition law forbids W — states with S > 1. However, we observe that the state |[gg] ' } ® I'&^V) i s 
a singlet and \&a F ) is a W = state. It has the right quantum numbers and, in principle, it could have a 
non-zero overlap with |\I>o(iV, U)). Nevertheless, the matrix elements of W between this state and the ones in 
the Ansatz are zero. This is the reason why we have dropped it in the expansion (43) . 
The three states in Eq.(43) are eigenstates of the kinetic-energy operator Ho: 

Ho\g°) ® |<0 = -2%°> ® lO 



Ho\g°) ® |<^° > = 2t\g°) ® |<^° ) 

1 ' m >®|*^ m > =0 



i 

(-) 



m— — 1 



(44) 



Expanding c^jC^f occurring in W in a, 6, 5, g operators of Eq.(7), we may take advantage from the fact 
that \$°/ F ) and m = 0,±1, are W = states. Indeed, since one cannot annihilate g or g over them, 

taking for instance x € A, 



(45) 



then, multiplying by <?j<?j, and the like, one obtains that the contribution proportional to 

^2 ^^Rj^iaj^oiN, U)) 



(46) 



yield nothing. Of course, if x £ £>, by the same reasoning we get 

^2Ri,a,Rj,xbiib jr \^ {N,U)) = 0. 



(47) 



Hence, we can write 



W = W A + W B 



with 



and 



W B :=U 



x=N+l 



2N 



2^ (SI + Sl)(ST +5t) + 



(-Si +fl , l)(-ST + St) + 




1 

] #i,*[(Sl + SiHt + a *l(ST + St)] 



(48) 



(49) 



1 N ~ 1 



(-ST + 3l) 



(50) 



where := means that the two sides are equivalent when acting on |^o(^V, U)). Next, we transform the remaining 
two operators, using Eq.(8). We note that the terms containing the sequence g^g^g a and similar terms 
with b instead of a and/or some g replaced by g vanish by symmetry; this is evident since X^=i Ri > x = f° r 
any i = 1, . . . , N — 1. Similarly, the sequences like a^g^gg also do not contribute. The sequences like a^a^gg 
are diagonal in the orbital indices of the two a 1 " operators, since y~) N _ 1 Ri, x Rj, x = 8ij\ hence, they annihilate 
|^ f o(^V, U)) ( recall that in this state all the a orbitals are singly occupied, as we can see by direct inspection 
of the expressions for |$^) and \&af), m = 0, ±1, of Section IV). By a particle-hole transformation one 
can show that the remaining terms that create pairs in the shell Shf, like a^a^ga, etc., have no effect, as in 
Eqs. (46-47); we postpone the proof of this last spectacular cancellation until Appendix A. To sum up, no term 
in W can change the number of particles in Shf, that is, that number is independent of t and U\ This proves 
the remarkable property that we have named Shell Population Rigidity: it holds for any finite TV and it is not 
specific of the thermodynamic limit. We shall see below that the Shell Population Rigidity characterizes the 
ground state and a suitable subspace of the total Hilbert space. In Fig. 2 we have drawn a physical picture 
of the cancellation among diagrams which do not preserve the number n a + n b . We can say that the W = 
states entering in our Ansatz are stable with respect to the Hubbard interaction W also in the presence of 
other particles in the system. This fact allows to exactly solve the Schrodinger equation. 
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(g, g) (g, g) 

w 



(a, b) 



W 



w 



(a, b) (g, g) 



(a, b) 



W 



= 



(8.8) 



(a, b) (a, b) 



(a,b) 



(8.8) 



(g.g) 



FIG. 2. Cancellation among the scattering amplitudes in the ansatz of Equation (43) which do not preserve the 
number of particles in the shell Sh f . 



Taking these cancellations into account one obtains 

w A ~ ^(g\+g\)(9l + g\)(9i + 9i)(9T+9x) + 



AN 

U 
2N 



(51) 



and 



U 



w E := j^(-g\ + g\)(-g\ + g\)(-9i + §i)(-9i + Sr) + 



AN 



+ ^fE [ & lT(-5l+P i t ) + (-5 T t +ff T t )fcl i ] [(-ffl + 5l)&<T+6ii(-ffT+ffT)] 



(52) 



With W written in the transformed picture one can calculate W\^o(N, U)). After very long, but standard 
algebra one finds 

W\g°) ® |<0 = u {^V> ® |<0 + 2^°> ® |<0+ 

JV-1 

i=l 
1 JV_1 

+ [sjffl ("a - n u) + g\g\(-nh + ntt)] \$a°f) + 

i=l 

N-l ~) 



and 



JV-l 

+ ^ E l^] 1 ' 1 ) ® ^ Sr")|<0 + 

i=l 
1 JV_1 

+ 2iV E [ffIfl{(~"iT + n ?T) +3{ffl(^a ""a)] \^a°f) + 

i=l 

JV-1 ~) 

+ ^Ei^] 1 '~ 1 >®(-^ +a +Oi<o ■ 
i=i j 

while for the singlet E^-iHHkw] 1 '" 1 ) ® I^af™) onc obtains 
i r i 



(53) 



(54) 



E 



(55) 
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By means of the identities in Eqs. (38-39-40) one can write Eqs. (53-54) in a compact form 
f 2N - 1 , 



W\g°) ® |<0 = [/ ■ 



2iV 



and 



W\f) ® |<0 = t/ j^V) ® l<0 + 2^15°) ® l<0 - ^ E (-) m |[33l 1 ' m ) ® l^ m )| 



(56) 



(57) 



Let us now consider the second row in Eq.(55). If our Ansatz is correct, the state in the square brackets 
must be proportional to l^ ^}. In this way one could close the equations and find an exact eigenstate. In 
Appendix B we prove that this is the case and in particular that 

JV-l 

E 

i=l 

Hence, Eq.(55) yields 
i 

w e (-r\\39] 1 n»\^ m ) = u 



(S+ a - S+ 6 )|$^- 1 ) + J=« T - nl - < + n^lO + (-S- a + S- b )\<f>Y F ) 



= V2(N 2 - 1)|^° > . 

(58) 



N + 1 



N 



E (-ri[^] 1 ' m >®i^ m >+ 



+ n/2 



m= — 1 

r2 



TV 



2iV 



V>-i<?°»®io 



(59) 



This result, together with Eqs. (56-57) and Eqs. (44), allows us to reduce the Schrodinger equation (^Hubbard — 
E)\' l ii i )(N, U)) = to the diagonalization of the matrix 



H{N,U) 



-2t + 



-U 



u 

2N 



^ 2(N . 2 ~ 1) u \ 



u 

2JV 



2t + 



2N-l r r ^2(N2-1) 

O AT U O AT 



U 



(60) 



N + l TT I 

\ 2Af 2JV N ' 



where the relations 



and 



[<5°l ® <0] [lff°> ® l<0] = [<5°l ® <0] [ls°> ® l<0] 



(61) 



(9°\®{* a° f \w E H m l[33] 1 ' m )®l$^ m )= E (-ri^n^ilggfnwig )®^), 

m— — 1 m— — 1 

1 1 

{g°\ ® ^ (-) m j[ 5 3] 1 ' m ) ® fe m > = E (-) m < $ AF m | ® <[ss] 1,m |W|s°> ® |<^° > (62) 

m— — 1 m— — 1 

have been used to express =-i (*af I^af ) m terms of { < l >( J i^ | 'I'a^) ■ 

The thermodynamic limit iV — * oo is well defined and non-trivial if we rescale t in such a way that Nt = 
t =const (which corresponds to the case of a fixed energy gap between the first two energy levels). The 
Hamiltonian matrix in Eq.(60) becomes 



with eigenvalues 



-2t + U U/V2 
H(oo,U)=\ 2t + U -U/V2 
U/V2 -U/V2 U 



£o = U, £± = U ± A, A = ^U 2 + At 2 



(63) 



(64) 
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The (unnormalized) eigenvector corresponding to the lowest eigenvalue £- gives 



|*o(oo,£0) = 



(j-^(2+j))b°)+ji/>]®i^>+V2(2+f) X)(-niff5] 1,m >®i*^ m >- , 

m=-l (65) 



So far we have proved that the Ansatz in Eq.(43) gives three exact eigenstates of the CBG-Hubbard model 
at half filling. The eigenvalue of lowest energy £- reduces to the non-interacting ground state energy for 
U ^ 0, while the associated eigenstate (65) to the state \g°) <g) |$^). On the other hand, when U — » oo the 
original Hamiltonian can be mapped onto what we can define the CBG-Heisenberg Hamiltonian: 



^Heisenberg = — J^( S * ' S V ~ j) • ( 66 ) 



U 4 



This model is exactly solvable and the ground state is obtained by projecting onto the singlet the state where 
N particles lie on the ,4-sites with spin up and the remaining JV on the £>-sites with spin down (Neel state): 
one finds 

N 

|* (iV, co)) = ]TR P [J al P(x+N) |0> (67) 

P x=l 

with aX t y — "^75( c i| c y| — c lj. c yt)' The ground state energy is —2t 2 /U for AT — » oo and is equal to the first order 
approximation of £_ in the small parameter (t/U) 2 (the same conclusion holds for any finite AT, although it 
is more tedious to prove). Furthermore, by direct inspection of the lowest energy eigenvector of H(N, U) one 
can show that |^o(iV, U)) becomes the state in Eq.(67) when U — > oo. Therefore, |^'o(A r , U)) reduces to the 
ground state of the CBG-Hubbard model in the two opposite limits U — > and U — > oo. 

To prove that the lowest energy eigenstate of H(N, U) is the unique ground state we are looking for, one can 
exploit the ground state uniqueness of the Heisenberg Hamiltonian proved in Ref. [26]. Since |^'o(A r , oo)) is the 
ground state of the Heisenberg model and since the ground state of the half-filled Hubbard model is unique, 
a level crossing for some value of U would be required if \^o{N,U)) were an excited state, contradicting the 
uniqueness. 

In conclusion we have proved that the Ansatz in Eq.(43) with (75,75,70) given by the lowest energy eigen- 
vector of the matrix in Eq.(60) is the exact ground state of the half- filled Hubbard model defined on the CBG 
of size |A| = 2 AT and repulsion U. In the next Section we shall calculate some physical quantities, as the 
spin-spin correlation function, and we shall prove that the particles are antiferromagnetically ordered. 



VI. RESULTS AND DISCUSSION 

The exact ground state solution for arbitrary but finite A^ allows to study the ground state energy Eq(N, U) 
as a function of A^ and U. Taking tN = t = 1, in Fig. (3) we have plotted Eo(N, U) in the range N = 1, .., 10 
and < U < 20. 




FIG. 3. Ground state energy Eq(N, U) of the half-filled Hubbard model on the CBG as a function of the on-site 
repulsion parameter U and the number of lattice sites N in each sublattice. The hopping parameter has been chosen 
to be t = 1. 
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The figure shows that Eo(N, U) is a monotonically increasing function of N and U, that is Eo(N + 1, U) — 
E (N) > and dE (N, U)/dU > 0. The limit N -> oo at fixed U is given by £L in Eq.(64): 



lim lim £? W (|A|,t,t/) = (7 - \JU 2 + 4t 2 < (68) 

| A| — >oo (3 — too 

where E^'(\A\, t, U) is the internal energy of the system at inverse Temperature (3 and number of lattice sites 
|A| = 2JV. On the other hand, one can calculate the internal energy at zero Temperature in thermodynamic 
limit: E (I3 \\K\ -> <x,t,U) = Qp + k B S//3 + fj,\A\n with Qp from Eq.(16) and S the entropy. Following the 
procedure of van Dongen and Vollhardt one gets 

.. r-rffli/iAi , m I — 2i = — 2Nt < if t is not rescaled ,„„. 

jELj&a El ^ m ' i ' D) = \o<i) iff~i/|A| ■ (69) 

In the case t=const the two limits commute and the trivial thermodynamic behavior of the model originates 
from the infinite gap between the lowest level and the zero-energy levels: the two particles in the g orbitals are 
completely decoupled from the dynamics of the system. On the other hand, if t is rescaled the gap remains 
frozen and taking |A| — > oo first, we can only say that the ground state energy is O(l), but we cannot predict 
the exact amount within the scheme proposed by van Dongen and Vollhardt. This is a consequence of the 
fact that in the limit |A| — > oo, we retain only the extensive contributions to the internal energy. 

The ground state energy Eq.(68) is always higher than that of the non-interacting case and hence non- 
trivial correlations survive in the thermodynamic limit. To understand the nature of these correlations we 
have calculated the ground state average of the number of doubly occupied sites: 

— W(N U) 1 1 

D(N, U) = v ' ' = -(* (JV,l7)|W|*o(JV, U)) = jj{Vo{N, *7)|tf H ubb ar d - H \^o(N, U)) = 

= V E (N, U) + -[ 7s 2 (iV, U) - j 2 s (N, U)] (70) 

where 7 S and 7g are the first two components of the normalized ground state vector, see Eq.(43). In Fig. (4a) 
we have plotted the trend of D(N,U) as Af increases for different values of U. As expected D(N,U) is a 
monotonically decreasing function of U, 8D(N,U)/dU < 0, approaching zero for U — > oo [where the exact 
ground state reduces to the antiferromagnetic Neel state, see Eq.(67)]. Nevertheless, D(N,U) shows a non- 
trivial behaviour as N increases at fixed U values. In the weak coupling regime, U <fC the number of doubly 
occupied sites grows as N becomes larger and larger converging to a finite value when N — > oo. The opposite 
trend is observed in the strong coupling regime, U 3> t, where D(N, U) decreases as N increases. In the 
intermediate regime U ~ t, D(N, U) is an increasing function of N for small N, but becomes an decreasing 
function for large N. Hence, there is a critical value U C (N) where the analytic continuation of D(N, U) to 
real N verifies 



\8UdN J 



= (71) 

U C {N) 



D W 




FIG. 4. a) Trend of the ground state average of the number of doubly occupied sites versus N in the range 
1 < N < 20 for different values of the Hubbard interaction parameter U =0.1, 1, 2, 3, 5, 10. b) Ground state average of 
the interaction term W as a function of U in the range < U < 10 for three different number of sublattice sites, N = 
1, 2, 10. The hopping parameter has been chosen to be t = 1 in both cases. 
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In Fig. (4b) we report the ground state average of the interacting term W as a function of U for three 
different values of N. There is always a critical value for U where the repulsion is maximized. Moreover, 
limjv^oo W(N, U) = Woo (U) 7^ for any finite U and the system cannot avoid double occupation neither in 
the thermodynamic limit. In the limit U — * oo, W(N, U) — » according to the Gutzwiller Ansatz [27]. 

Next, we have calculated the spin-spin correlation function 



G spill (x,y) = (t> (N,U)\S z x S z y \t> (N,U)) = lim G spin (x,y; f3, h = 0, |A|) = lim lim G Bpin (x,y; p, h, |A|) 



(72) 



where G Bp in(x,y; (3, h, |A|) was defined in Section III, S£ = \{n x i — n x i) and |*o) is the normalized ground 



state 



l*o> 



+ ±^y, (-r\m Um ) ® i*i; m ) , (73) 



m— — 1 

with A/? = (*^f|*^>, A/? = <$^™ |*^>, m = 0,±1, and 7s + t| + 7o = 1- Due to the Sn ® Sn ® Z 2 
symmetry, G sp i n (x,y) can be written as 



{Go x = y 
G on x € A (x G B) and y G A (y G £?) 
Goff a; € .4 (a; G B) and 1/66(16^) 



(74) 



and by exploiting the sum rule gA G Bp i n (x, y ) = Go + (N — l)G on + NG g = 0, one can express G on in 
terms of G Q ff and Go: 



N r Go 

"(-Toff 



iV — 1 N — 1 

The problem is then reduced to the evaluation of Go and G «- We have 

G„ = Gsp^x) = <*o|(^) 2 |*o) = ± ^<*o|(S|) 2 |* ) = 



(75) 



— ^(* |n xT - 2n xT ri a;i |*o> = -(1 - -^) 



(76) 



16A 



More effort is needed to calculate G g. Expanding the c operators in g, g, a and b operators of Eq.(7) and 
taking into account that (^2 x£A S x ) 2 \^fo) = 0, we get 



(77) 



with (— ) CT = +, — for a =|, J.. After long but standard algebra one obtains 



-A^[(^-7 5 ) 2 +27o 2 ]-Vf 



A/oM 



3A? 



and exploiting the symmetry <ij — ► bj and bj — > a,, which implies |$ A ' F s ) — ► (~) S | ( 3? y iV S )> 



(78) 

(79) 
(80) 
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M 2 



(81) 



where in Eq.(79) we have used the identity Eq.(39). Finally we recall that the W = states \^^ F s ) arc 
eigenstates of the square of the total spin operators of each sublattice S\ = (S za ) 2 + ^(S +a S~ a + S~ a S +a ) 
and Si ee (S zb ) 2 + i,(S +b S- b + S- b S +b ) with eigenvalue (2=1) (£1=1 + 1) = i (N 2 - 1); hence 

s 



(*AF S \(S Za ) 2 \<S> S AF S ) 

Mg = -S 



N 



12 



-a 



Substitution of Eqs. (79-80-81-82) in Eq.(78) yields 



iV 2 



12 



1 . 1 

8 



+ ^Klg - is) 2 + 2 7o] - tj7o 



N 2 



18 



-7o(79 -7s) 



(82) 



(83) 



where we have taken into account Eq.(62) to evaluate the ratio JV1/A/0 = y (A^ 2 — l)/3. 




u 




FIG. 5. a) G n versus £/ in the range < £/ < 20 for three different values of the number of sites N — 2, 3, 10. 
b) G a ff versus U in the range < U < 20 for three different values of the number of sites N = 2, 3, 10. The hopping 
parameter has been chosen to be t = 1 in both cases. 



In Fig. (5) we report the trend of G on and G ff versus U for three different values of N = 2, 3, 10. According 
to the Shen-Qiu-Tian theorem [28], G on is always larger than zero while G a g is always negative. Now consider 
the ground state average of the square of the staggered magnetization operator 

m z AF = ±.(i S f \[^2e(x)S z x ] 2 \^ ), e(x) = 1, -1 for x e A, B . (84) 

The Shen-Qiu-Tian theorem implies that each term in the expansion of Eq.(84) is non-negative. We emphasize, 
however, that for \A\ = \B\ this does not imply that m 2 AF is an extensive quantity! Remarkably, our results 
show that m 2 AF — Go + (N — l)G on — NG a a is extensive for any value of the on-site repulsion parameter U 
and provide an example of antiferromagnetism in the ground state of an itinerant-electrons model of arbitrary 
size. 

We also observe that the thermodynamic limit yields a non- vanishing result: 

lim Go = 7, lim G on = - lim G off = --, (85) 

A|— >oo 4 |A|— >oo |A|^oo 12 

or 

lim lim lim \G ap i n (x, y; j3, h, |A|)| = -7- for x ^ y . (86) 

|A|— >oo /3^oo h— >0 12 

By comparing Eq.(86) with Eq.(23), we see that the lim^^o lim^^o do not commute with the IimiA|->oo- 
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VII. SUMMARY AND CONCLUSIONS 



We have introduced the CBG-Hubbard model and explicitly written the ground state wave function at half 
filling for any repulsion parameter U. The topology of the graph A allows (2N — 2)-body eigenstates \<&af S } 
of the Hamiltonian which are free of double occupation. They can be obtained projecting the determinantal 
state |<&af) on a given spin-S 1 subspace. |<&af) has the antiferromagnetic property, that is the map A <-> B is 
equivalent to a spin-flip. 

We have shown that the ground state has a fixed number of particles in the zero-energy one-body shell Shf, 
independent of t and U. This very remarkable Shell Population Rigidity holds in any finite-size system and 
it provides the closure of the equations. This implies extra-conservation laws in a suitable subspace of the 
full Hilbert space including the ground state. Only the singlet |$^} and the triplet |<I>^™), m = 0, ±1, are 
involved in the ground state expansion. Qualitatively, we may say that according to Eq.(65) the particles in 
the shell Shf manage to avoid the double occupation. 

We have calculated the spin-spin correlation function and shown that the ground state exhibits an antifer- 
romagnetic order for any non-zero U even in the thermodynamic limit. The kinetic term induces non-trivial 
correlations among the particles and an antiparallel spin configuration in the two sublattices turns out to 
be energetically favoured. Therefore, the model contains the basic ingredients to understand how the subtle 
competition between the derealization induced by Ho and the localization induced by W may give rise to a 
magnetically ordered state even outside the strong coupling regime (where the Hubbard model is equivalent 
to the Heisenberg model). On the other hand, within the scheme proposed by van Dongen and Vollhardt, 
the model is well described by a paramagnetic Hamiltonian; the difference stems from taking the thermody- 
namic limit before the limit of zero Temperature, because then the kinetic term completly decouples from the 
interaction and loses any role in the physical response functions. 

The present formalism lends itself to solve more realistic Hubbard Hamiltonians with an increasing number 
of negative and positive energy levels. Higher-spin projections \^^p S ) are involved in the ground state 
expansions; the results will be published elsewhere. We are currently investigating group-theory aspects of 
this model and its extensions. 

Finally, we recall [23] [24] [25] that the standard Hubbard model on a N x N square lattice and periodic 
boundary conditions also has |$ A '^5 s )-like eigenstates. It could be that the Shell Population Rigidity holds in 
this case too, but the proof is lacking. 



APPENDIX A: SHELL POPULATION RIGIDITY 

We have seen in Section V that the most part of the contributions which do not preserve the number of 
particles in the shell Sh / yield nothing. To prove the Shell Popoulation Rigidity property we need to show that 
the sequences like a 1 a) ga and similar terms with b instead of a and/or g replaced by g vanish or in formula 

J2 Yl R i,x R i\*^,*4Al a ™\®AF S ) = ( Al ) 

x£*A ijm 

and the like with a replaced by b. This last spectacular cancellation can be seen by particle-hole transforming 
Eq.(Al). Under a particle-hole transformation c x<r — > c\ a and the W = states \&Ap S ) ~ ► \g°g°) <8> |$!i? Ms )i 
modulo a phase factor. Hence, Eq.(Al) is equivalent to 

^^/.', , /.', /.',., <I'h U - = . (A2) 

x£A ijm 

Let the spin of the operator a) ma be up for the sake of definiteness (the same reasoning holds in the case of 
down spin). The l.h.s. of Eq.(A2) can be rewritten as 



xGA- ijm j |_ x i 



ajl\* S A? s ) (A3) 



where we have exploited the lack of double occupation of the W = state, see Eq.(45). Therefore, we are left 
to prove that 

JV JV-l 
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We recall that R is an (JV— 1) x iV rectangular matrix whose rows are orthonormal vectors which are orthogonal 
to the A''— dimensional vector (1, 1, . . . , 1, 1). It is a simple exercise to verify that 



( 75 -75 ... ^ 

... 



R = 



1 1 2_ 

\/6 \/6 s/3 



x > m + 1 



\ ^/AT(JV+1) ^/]V(iV+l) ^JV(iV+l) ^/iV(]V+l) " ' V J 

is a correct choice and that Eq.(A4) is identically verified. 



Rm,* = { -v^T x = m + l 

x<m (A5) 



APPENDIX B: PROOF OF EQ.(58) 

To prove Eq.(58) we shall use the definitions (29-35-36) for the triplet W = state |$^™), m = 0, ±1. The 
first term yields 

- sjf*)!*^ 1 ) = E(-) p E(^ +a - E II i°) - 

k p k j i^j 

y+a _ Q +b W (-l)t TT t 



= E(-) p E(^ -^))^n-U)i°) + 
+ E(-) p E« t )E^ +a -^)K^) II *Ui°> 
= ^Eh p E^o)I1-U)I ) + 

P 3 i& 

- ^B-)'EE*^) n *Ui°> = 



r (-l) t T (+ 1 ) t 

while second term yields 

E^T -^fei - +«t l )i , i > AP> = E^) p E^t - A ki - npwi + A p(m)^2 T j PU)'[[' 7 tp(i)\ ) = 

k p k j 

= E ( ~ )P E ( ^T - ™i! - fi b P(m + np 0)i )T^/ a) <rt p(j) |0) + 

P 3 i& 

+ E^ P E T ™0)E^T - ft fcl -»P(fc)T + A P(fc)i)°"l,P(fc) II = 
P j fc^i <^j,* 

= 2 E(-) p E-U)Il-U)i ) + 

+ 2 EH p EE^US( fe )n-U!o)- 

P j fe^j i^i,k 

-2 ( iv- 1 )i^) + 2E(-) p EE^%)^, II • ( B2 ) 

The third term can be computed following the same steps of Eq.(BI) and the final result is 

E(-*r + s- k b )\^ F ) = V2( N - wr F ) ^Eh p E E ^Ste'S, II <^i°> • 

k P j k^j i^j,k (B3) 

Denoting by [..] the left hand side of Eq.(58), from Eqs.(Bl-B2-B3) one obtains 

[..] = sv-2(n - 1)1*50 > - ^E(-) P EE Hr^^S) - ^U°4> + #2W:$>] II *Ui°>- 

P j k^j i^3,k 

(B4) 
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Let us consider the first and the third term in the square bracket of the right hand side of Eq.(B4). We have 
Vf-l P ' M~ 1)f r (+1)t +r (+1)+ r ( - 1)f 1 TT ^ 10) - 



i^3,k 



= -£(- 



r (-D T r ( +1 ) T +r (+ 
T j,P(k) T k,P(j) ^ 'j,p 



■(k) T k,p(j)] n a ip(«)i o ^ 



i^j,k 



where the permutations P and P' satisfy 

P'(fc) = P(j), P'(j) = P(fc), P'W = P(i) Vz j£ j, fc ; 
hence (— ) p = — (— ) p . Substituting Eq.(B5) into Eq.(B4) and taking into account that 



,(0)+ fo)t 



r (+i)t (-i) f 



4 



' j,P(k)'k,P(j) ^ 3,PU)'k,P(k) T 'j,P(k)'k,P(j) - u j,P(j) u k,P(k) 



one obtains 



[..] = 3V2(N- l)|cO + V2^(-) P ^^4 P( , } < P(fc) J] at pw |0) = ^(iV 2 - 1)|<0 

P j k^j i^3,k 



that is, Eq.(58). 



(B5) 



(B6) 



(B7) 



(B8) 
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